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Two atomic quantum dots interacting via coupling to BECs 
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Abstract. - We consider a system of three weakly coupled Bose- Einstein condensates and two 
atomic quantum dots embedded in the barriers between the condensates. Each dot is coupled to 
two neighboring condensates by optical transitions and can be described as a two-state system, 
or a pseudospin 1/2. Although there is no direct coupling between the dots, an effective interac- 
tion between the pseudospins is induced due to their coupling to the condensate reservoirs. We 
investigate this effective interaction, depending on the strengths of the dot-condensate coupling 
T and the direct coupling J between the condensates. In particular, we show that an initially 
ferromagnetic arrangement of the two pseudospins stays intact even for large T/J. However, an- 
tiferromagnetically aligned spins undergo peculiar "breathing" modes for weak coupling T/J < 1, 
while for strong coupling the behaviour of the spins becomes uncorrelated. 



Introduction. — Atomic quantum dots (AQDs) are 
a relatively new topic, referred to as nanobosonics in the 
context of cold atoms. An AQD constitutes a lowest 
atomic level of a tight trapping potential which can be 
maximally singly occupied due to a large repulsive inter- 
action energy between the atoms. In analogy to nano- 
electronics the question arises whether one can couple an 
atomic dot to large bosonic reservoirs, i.e. "leads" and in- 
vestigate particle transport properties in such a combined 
system. Another interesting problem is the control and 
manipulation of the single spin behaviour, since such a 
quantum dot can be described in terms of a two-state sys- 
tem (the occupied state is equivalent to a "spin-up" , and 
' the empty state corresponds to a "spin-down" ) . 

An interesting variant of the optical coupling between 
an AQD and a uniform superfluid reservoir was consid- 
ered by Recati et al [1]. Provided the dot atoms and the 
atoms of the superfluid reservoir are of different hyperfinc 
species, one can (by an external laser) induce Raman tran- 
sitions between dot atoms and those of the condensate. 
At low energies the spin degrees of freedom couple to the 
lowest excitation modes of a superfluid bath - phonons, 
and therefore the combined system can be mapped onto 
a spin-boson model, and the dissipative behaviour of the 
spin can be explored [1]. From quantum optical point of 
view this kind of Raman coupling and its consequences 
were numerically investigated by Zippilli and Morigi [2]. 




Fig. 1: System of two atomic quantum dots coupled to three 
Bose-Einstein condensates in a three- well potential Vbec- In- 
side the barriers between the wells atomic quantum dots (single 
atoms of a different from condensate bosons hyperfine species) 
are embedded. T is the optical coupling between the dots 
and condensate, and J is the Josephson tunneling between the 
wells. 



In previous works by Bausmerth et al. [3] and Fischer et 
al. [4], a single AQD coherently coupled by optical tran- 
sitions to two BEC-reservoirs with finite number of par- 
ticles has been considered. A direct Josephson tunneling 
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Fig. 2: Rotation frequencies (Eqs. (10) and (11)) of the pseu- 
dospins for symmetric initial conditions in the case of weak 
coupling. We took iVi(0) = N 3 (Q) = 500, N 2 (0) = 1000, 
U = and the relative coupling of the dot to the condensate is 
T/J = 10. The frequencies uj 3 x and ui 3 y are expressed in units of 
2J. It is clear that in both FM and AFM cases the spins nu- 
tate. Plot (b) shows a zoom in details of plot (a) on a scale, on 
which the differences between the three curves in plot (a) can 
be distinguished. Red and black curves show the frequencies of 
the two initially antiferromagnetically aligned spins, while the 
green line is for initially ferromagnetically aligned spins. The 
detuning energy is fixed to H5/2J = 1. 



between the two BECs has also been taken into account. 
One should note here, that a Bose Josephson junction ex- 
hibits an unusually rich dynamical behaviour [5-7] , which 
was also confirmed experimentally [8]. 

In previous works [3] and [4] it was observed, that one 
should distinguish between the two important limiting 
cases: (i) the weak coupling regime, in which the cou- 
pling between dot and condensate T is smaller or of the 
same order as the direct Josephson coupling J between 
two condensates, and (ii) the strong coupling regime, in 
which the coupling between dot and condensate is domi- 
nating, so that the tunneling between the two condensates 
occurs predominantly through the dot channel [4] . 

Since in the weak coupling regime the effect of the AQD 
on the tunneling between the condensates is very small, 
it is the favorable regime in order to investigate the be- 
haviour of the pseudospin under the influence of an effec- 
tive time-dependent field, induced by the condensates dy- 
namics. It turns out, that the pseudospin behaves in this 
case as a quantum top and undergoes multi-frequency ro- 
tations depending on the initial conditions and the mean- 
field interaction parameter [3]. 

In the strong coupling limit (T> J), it was shown that 
an atomic quantum dot, surprisingly, can cause large par- 
ticle imbalance oscillations between the wells similar to the 
Josephson effect, which are characterized by the two fre- 
quencies. For non-interacting condensates these frequen- 
cies could be derived analytically [4]. 

In this work we discuss the three-well condensate with 
two AQDs embedded between the wells (fig. 1), coupled to 
the condensates by optical Raman transitions. We show 
that, although the AQDs do not directly interact with 
each other, their coupling to the condensates induces an 
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Fig. 3: Relative orientiaton sis*2 = cos(7) versus time for dif- 
ferent coupling constant T/J. The initial conditions are the 
same as in Fig. 2 and initial alignment of the two pseudospins 
is AFM. One can observe the "breathing" modes, whose period 
decreases for growing T/J. 



effective coupling between the two pseudospins, associated 
with the two dots, which can be studied numerically. We 
show below that symmetric initial conditions favour the 
ferromagnetic (FM) coupling between the dots, while an- 
tiferromagnetic (AFM) spin alignment does not survive 
the increasing coupling T/J. These features are especially 
clearly observed in the strong coupling limit. 

Model. — We investigate the setup depicted in Fig. 1. 
The condensate in a three-well potential V BEC is described 
within a three mode approximation, which is just a gen- 
eralization of the two-mode case [5,6], so that the wave- 
function is ^(r, t) = Y^i=i 4>i{ v )^ i where the </>j(r)'s 
are the normalized local mode solutions for well i, and 



*i(«) = yjVj(t)e iS<(t) . 



(1) 



Here N% is the number of particles in the i-th condensate, 
and 0i is its phase. 

The Hamiltonian of our system reads 



H = ^^|*i| 2 + ^^|*i| 4 -J(^*a + *S*i) 

i=l i=l 



(2) 
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{ai j) ± ia^>)/2, where ag' , atf>, and a 
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where a 

are the Pauli matrices describing the pseudospin degrees 
of freedom of the first (j = 1) or the second (j = 2) atomic 
quantum dot. Each quantum dot is described by a state- 
vector = (^ j d \a^ \^/ 3 d ) with j = 1,2 and ^ J d being a 
two-state wave function of the j-th dot. This formalism is 
applicable in the limit of large interaction on the atomic 
dots. 
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Fig. 4: Behaviour of the pseudospins for asymmetric initial 
condensate sizes in the case of weak coupling: iVi(O) = 10000, 
7V 2 (0) = 1000, AT 3 (0) = 500, interaction is fixed to UN /2J = 
1, and the relative coupling of the dot to the condensate is 
T/J = 0.1. Plot (a) displays the frequency evolution in the 
Lo x -LOy-plane, while in plot (b) the behaviour of the pseudospins 
is shown. The spins were initially FM aligned. The black curve 
denotes the left pseudospin, and red curves correspond to the 
right pseudospin. 



For simplicity we consider a symmetric system: the 
Josephson couplings between the condensates are the 
equal: J± 2 — J23 = J, defined through 



(V^(r)V0 i+1 (r)) 

2m 



+&(r)V r MO (r)0 i+1 (r) , i = l,2 



(3) 



All four dot-condensate couplings are expressed by a 
single parameter T 

T = hn R f dr<t> j (r)4(r) = hn R f dr^ i+ i(r)^(r), 



(4) 

where j — 1,2, and the spatial wave function of each dot 
<^(r) is normalized to unity, flu is the Rabi frequency of 
the Raman transition. The Rabi frequency is an externally 
controllable parameter, that can be tuned to any desir- 
able value. Spontaneous emission is suppressed by a large 
detuning from the excited electronic states, which is ab- 
sorbed into the effective dot energy HS [1] . Another impor- 
tant parameter in the problem is the standard mean-field 
two-particle interaction between the condensate atoms 



U = g / dr|^(r)| 4 , 



(5) 



where g = 4ira s /m with a s being the s— wave scattering 
rate. The zero-point energies are taken equal to each other 
E i = I [-£ I V <M r ) I 2 +\Mr)\ 2 V BEC (r)\ dr = E, and 
in the actual calculation we assume E = 0. 

We can now derive equations of motion for the conden- 
sates wave functions, which take the following form 



id t ^>2 = {E + U\^ 2 \ 2 )^2~ J(*l + * 3 ) 

+ T( S « + S ? ) ), 
id t ^> 3 = {E + U\V 3 \ 2 )V 3 - JV 2 + Ts {2) 



(7) 
(8) 



The equations of motion for the pseudospin components 
of the two quantum dots are expressed as two Bloch equa- 
tions 

hdtSV) = x ^ . (9) 

The Bloch equations for the two spins are coupled through 
the rotation frequencies 



and 



t« = 



2T(m! + m 2 ) 

-2T(3*i + 3^2) 
-US 

2T(m 2 + 3^3) 
-ns 



(10) 



(ii) 



which are time-dependent due to a the coupling to the 
condensate dynamics Eqs. (6)- (8). 

We solve this set of equations numerically for different 
initial conditions, and analyze the behaviour of the con- 
densates and the dots in both regimes of weak and strong 
coupling. In case of the condensate we are interested in 
the behaviour of normalized particle imbalances 



yvift)-yv 2 ft) 

N 



n 23 



N 2 (t)-N 3 (t) 
N 



(12) 



where N = iVi(0) + A^ 2 (0) + yV 3 (0). Note, that conserving 
quantities are N tot 
(s«) 2 



iVo + + l)/2 + («r + l)/2 and 



(2) 



(4 j) ) 2 + (4 J) ) 2 + (^' ) ) 2 for j 



1,2. 



idt*i 



(E + C/|*i| 2 )*i - J* 2 + Ts 



(i) 



(G) 



Weak coupling limit. — The weak coupling limit im- 
plies that T < J (we note that, since a quantum dot can 
be only singly occupied, T can be also order of magnitude 
larger than J, but we still refer to this realization as the 
weak coupling limit). In this limit the effect of the AQDs 
on the dynamics of the Josephson tunneling between the 
condensates is expected to be insignificant [3]. We con- 
centrate therefore on the pseudospin behaviour. 

It turns out, that symmetric initial conditions, which 
means that the two outer wells are equally populated 
yVi(0) = N 3 (0), favour nutations of the two pseudospins 
(Fig. 2(a)), which remain parallel if they were initially 
FM aligned, and antiparallel if they were AFM aligned 
at t = 0. This is true for a very weak coupling T <C J, 
when the effective interaction between the spins is negli- 
gibly small. 

However, with increasing coupling, a certain correlation 
between the two pseudospins develops, and their dynami- 
cal behaviour starts being dependent on their initial align- 
ment. These features are to be observed in Fig. 2(b), in 
which part of Fig. 2(a) is enlarged. The two FM spins re- 
main parallel and their frequencies follow the same (green) 
curve. However, initially antiparallel spins do not remain 
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Fig. 5: Example of large amplitude particle oscillations be- 
tween the wells ni2 and niz in case of strong coupling limit. 
Left column (a) is for initial conditions Ni(0) = A^(0) = 
-^3(0) = 1000, right column (b) is for initial conditions Ni(0) — 
N 3 (0) = 500, N 2 (0) = 1000. For all plots the effective interac- 
tion UN Q /T = 0.01 and the detuning S/T = 100. Upper plots 
are for initial FM arrangement of the two AQDs, and lower 
plots are for initially AFM dots. 

antiparallel with time. The red curve in fig. 2(b) corre- 
sponds to the left dot, while the black curve corresponds 
to the right dot. 

In order to understand better this behaviour we ana- 
lyze the time evolution of the relative angle 7 between the 
spins, defined through 

S1S2 = cos 7. (13) 

The time-dependence of cos 7 for different coupling con- 
stants T/J is shown in Fig. 3. The angle 7 oscillates be- 
tween 7r and ~ 0.75-7T. In the following we refer to these 
oscillations as "breathing" modes. The breathing modes 
become slower for smaller coupling T/J and eventually 
disappear for T <C J. We conclude that the breathing 
modes are indications of the induced correlation between 
the spins. 

One can also consider asymmetric initial conditions: 
Nt(0) ^ N 2 (0) ^ jv 3 (o) as, for instance, in Fig. 4. In this 
case the dynamics of the two AQDs becomes very compli- 
cated, although not chaotic, as the Fourier transform of 
(cos 7) does not produce a typical white noise signal. 

Strong coupling limit. — The peculiarities, which 
we found out in the dynamics of the initially AFM aligned 
spins in the weak coupling limit, are expected to become 
even more pronounced in the strong coupling regime, cor- 
responding to T > J. In this case, the relevant time scale 
is set by T. It was previously shown that the ratio UN /S 
plays an important role in the physical behaviour of the 
system [4]. Thus, the most interesting regime, when the 
AQD serves as a quantum "shuttle" of the atoms between 



the condensates, is achieved for S ;§> UNq. Naively one ex- 
pects, that the oscillations induced solely due to the cou- 
pling of the condensates to the dot, would be always very 
small, however, surprisingly large amplitude Josephson- 
type oscillations between the wells were obtained [4]. 

In our case we also observe this property of AQDs 
(Fig. 5). For initially parallel spins we observe symmet- 
ric oscillations of the particle imbalances nn and 77,23 (see 
Figs. (5a) and (5c)). These are similar to the oscillations 
in a triple well potential which are always symmetric for 
symmetric initial conditions. Note, that these oscillations 
are induced purely by the dots. 

The pseudospins remain parallel, independent of the 
coupling strength between them, so that cos7(£) = 
cos 7(0) = 1. The frequency uj (which is the same for 
left and right spin) exhibits precession (Fig. 6(a)), which 
means that each of the pseudospins nutates. 

The situation is very different for initially antiparallel 
pseudospins. First of all, the dynamics of intcrwcll tunnel- 
ing is of different character and seems to be not very sen- 
sitive to initial conditions for A^(0) (Figs. (5b) and (5d)). 
Second of all, one clearly observes asymmetric features, 
which we associate with the uncorrelated spin behaviour. 
This can be also seen from the behaviour of uj^ and u)( 2 \ 
which is shown in Fig. (6b). For instance, one can no- 
tice, that the right (red) pseudospins oscillates faster. Fi- 
nally, the time-dependence of cos(7) is depicted in Fig. (7), 
where the uncorrelated dynamics of the two initially an- 
tiparallel (cos(7)(i = 0) = —1) pseudospins is presented. 
For small-amplitude oscillations the qualitative behaviour 
of the two AQDs is similar. 

Conclusions and Discussion. — We have discussed 
the problem of an induced interaction between two atomic 
quantum dots which are embedded between three trapped 
Bose-Einstein condensates. As it is a multi-parameter 
problem, it helps to distinguish two different regimes: a 
weak coupling limit (when the direct Josephson coupling 
between the condensates is dominating) and a strong cou- 
pling case, in which the tunneling between the condensates 
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Fig. 7: Time-dependence of sis*2 = cos(7) for the strong cou- 
pling case. Initial conditions are as in Fig. 5(b). 



occurs predominantly through the additional channels me- 
diated through the dots. 

Our main conclusion is that indeed we clearly observe 
the signs of an effective interaction, induced between the 
two dots due to their coupling to the same condensate sys- 
tem. For this interaction to appear, T should be at least of 
the order of J. Interestingly, this induced interaction pro- 
duces very different effects on the pseudospins, depending 
on whether they were initially FM or AFM aligned. 

It turns out, that initially parallel spins are stable with 
respect to coupling between the dots and the conden- 
sates. They remain parallel even in the strong coupling 
limit, while undergoing two-frequency behaviour, as we 
can judge from the precessions of Co. 

Initially antiparallel spins perform what we call "breath- 
ing" motion, which is characterized by small periodic de- 
viations of the relative spin angle 7 from tt. The breathing 
modes, however, do not survive the strong coupling limit, 
when the effective interaction between the two spins be- 
comes large. Instead, spins become rather uncorrelated. 

Our analysis has important consequences for the peri- 
odic model of the system under consideration, which is 
the subject of future research. First of all, in a periodic 
system one shall be able to locally address and modify the 
interaction between the neighboring spins just by chang- 
ing the occupation of each well. Second of all, by changing 
the coupling between AQDs and condensates (or optical 
lattice sites), one can control the different phase transi- 
tions in such a coupled system. For example, if an optical 
lattice was initially in the Mott state, a finite T should 
drive the system to a superfiuid state once some critical 
value is exceeded. At the same time spins of the AQDs 
will become correlated, and the question will be if this 
correlaion is associated with a phase transition. 

In the future it will be definitely worth studying quan- 
tum effects in such a model, as well as temperature ef- 
fects, which would lead to decoherence phenomena and 
temperature-dependent phase transitions. 
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